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■ Abstract 

(N 

This paper considers data rate limitations for mean square stabilization of uncertain discrete-time 
linear systems via finite data rate and lossy channels. For a plant having parametric uncertainties, a 
necessary condition and a sufficient condition are derived, represented by the data rate, the packet loss 
probability, uncertainty bounds on plant parameters, and the unstable eigenvalues of the plant. The 
results extend those existing in the area of networked control, and in particular, the condition is exact 
for the scalar plant case. 



*\ ! 1 Introduction 

Networked control systems have attracted much research interest in recent years [1,2]. In such systems, 
the communication between the plant and the controller is constrained due to the use of shared channels. 
Though modern communication channels may have large bandwidth, in the systems consisting of sensors, 
actuators, and other devices, each component can be allocated only a portion of it for the required real-time 
\ transmissions. Thus, it is important to study fundamental relations concerning the control performance and 

the communication constraints in such systems. 

One such relation is the limitation on the data rate for the stabilization of unstable linear systems. The 
seminal works of [3-5] have presented the minimum data rate for stabilization, and have shown that it only 
depends on the product of the unstable eigenvalues of the plant. Another fundamental relation can be found 
under the presence of packet losses. In practical channels, transmitted packets may be lost due to congestion 
or delay. In [6], it is shown that the maximum loss probability tolerable for mean square stabilization 
is described also by the product of the unstable eigenvalues of the plant. It is worth noticing that both 
limitations require more communication for plants with more unstable dynamics. On the other hand, the 
works [7, 8] focus on the two constraints of data rates and packet losses simultaneously and extend the results 
on limitations. It is interesting that while the work [7] studies the notion of stability with probability 1, 
the paper [8] is concerned with mean square stability. As a consequence, different bounds are obtained. 
Furthermore, various communication constraints are considered (e.g., delays, variable sampling periods, 
sharing the channel among multiple nodes), not only for linear systems but also for nonlinear systems (see, 
e.g., [9-11]). 

While in the existing works, it is commonly assumed that the exact plant model is known, when the 
plant is uncertain, fewer results are available. The work of [12] deals with a linear time-invariant plant with 
norm-bounded uncertainty; and in [13], a nonlinear scalar plant with stochastic uncertainty is considered. 
Both papers give sufficient conditions on the data rate for stabilization using a lossless channel. However, 
these limitations may contain some conservativeness, and it is still unclear at least how large the data rate 
should be when the plant is uncertain. 
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Yokohama, 226-8502, Japan. Emails: okano@sc.dis.titech.ac.jp, ishii@dis.titech.ac.jp. This work was supported 
in part by the Ministry of Education, Culture, Sports, Science and Technology, Japan, under Grant-in- Aid for Scientific Research 
Grant No. 23760385. 



1 



Plant 



Uk 



Vk 



Encoder 



7k- 1 



Controller 



Sk 



. Lossy 
\ channel; 



Decoder 



IkSk 



Figure 1: Networked control system 



The objective of this paper is to provide an answer to this question. We consider plants with parametric 
uncertainty. Though this class of uncertainty is different from those in the works mentioned above, it is widely 
adopted and various results are established (see [14] and references therein). To deal with the uncertainty, 
we follow the approach of [8], where the plants are given in the controllable canonical forms. Hence, the 
uncertainty can be considered as a natural extension of this approach. As to the communication constraints, 
following the setup in [8], we consider the channel is with data rate constraints and packet losses. 

We derive a necessary condition and a sufficient condition, which provide bounds on the required data 
rate and loss probability to stabilize the closed-loop system. The results have several features: First, the 
limitation given by the conditions is tight for the scalar plants case, which is not attained in [12,13]. Second, 
it generalizes the results in [8] in the sense that without uncertainty in the plants, the bounds coincide with 
that in [8]. Moreover, it is known that for plants without uncertainty, stabilization is possible with any data 
rate greater than the necessary bound by introducing a time-sharing protocol [4,5,8]. However, we show 
that this may be difficult in the uncertain case. 

The paper is organized as follows. In Section 2, we describe the setup of the networked control systems 
considered. Then, we present a necessary condition and a sufficient condition in Sections 3 and 4, respectively. 
The limitation for the system employing a time-sharing protocol is discussed in Section 5, and finally we 
provide concluding remarks in Section 6. 

2 Problem Formulation 

We consider the stabilization of a networked control system which includes a digital channel between the 
plant and the controller. At time k G Z + , the encoder observes the plant output y k € R, and quantizes it 
as an i?-bit signal s k € where R > 0. Here, the set S fl represents all possible outputs of the encoder, 
and contains 2 R symbols. The quantized signal s k is transmitted to the decoder via a lossy channel. Based 
on the received signal, the decoder generates the interval 34 C K, which is an estimate of the plant output 
y k . Finally, the controller receives this set 34 and provides the control input u k € K. Fig. 1 shows the 
connections among components of the system. 

Now we describe the details of each element in the system. The plant is an n-dimcnsional single-input 
single-output autoregressive system whose parameters have uncertainties and may be time varying: 

2/fc+i - ai,kVk - a 2 .kVk-\ a«,fe?A-n+i = u k- (1) 

Here, the initial value y of the output is bounded with a known interval \yo\ < Y , and yk = for k < 0. 
Each uncertain parameter a, ; fe is represented by the nominal value a* and the width > of the perturbation 

as 

a iik &Ai:=[a* - ei,a* + ei] for i = 1, 2, . . . , n, (2) 

where Ai represents the uncertainty of the ith parameter a^ k - 

The plant (1) can be rewritten in the controllable canonical form as 

x k+1 = A k x k + Bu k , y k = Cx k , (3) 
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where x k := [yk-n+i Vk-n+2 
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Let ^4* denote the nominal A-matrix, and let A^* be the product of the eigenvalues of A*: 

1 ••• 



A* :-- 











1 



, \ A *:=l[\i(A*) 



i=i 



where A^(-) represents an eigenvalue of a matrix. Assume that all eigenvalues of the matrix A k are unstable 
for the uncertain parameters in (2). In particular, this implies 



Kl - e n > 1- 



(4) 



The encoder Ek at time k quantizes the output yk € R to the i?-bit signal Sk € and then transmits 
it to the decoder via a lossy channel. Let the random variable jk represent the state of packet reception/loss 
at time k. The encoder knows the state of the previous step, 7fc-i, through the acknowledgement signal 
from the decoder. 

In order to deal with the uncertainty in the plant, we introduce some structure in the encoder. Specifically, 
it is realized by a finite level uniform quantizer. Such a class of encoders has been employed in, e.g., 
[5,8,12,13]. We will have more discussion later. Let 5at(-) denote the iV-level uniform quantizer with 
N = log 2 R and the input range [—1/2, 1/2]. It partitions the range [—1/2, 1/2] into N intervals of widths 
l/N as 



Qn{v) ■= 



if -I 



N- 



£<2/<-5 + Tr> ie{0, l,...,N-2}, 

tf |-£<y<§. 



Let the set of its outputs be Y, R := {0, 1, . . . , N — 1}. Then the encoder Ek is given by 



s k =E k (ytlo X ) = 



where j/q is the sequence y$ := {yo, . . . , yk} and similarly 7q _1 := {7o ; • ■ • j 7fe-i}- We use similar notations 
in the sequel. In addition, o k > is the scaling parameter of the input range. The above equation implies 
that the input range of the encoder is [— Ofc/2, crfc/2]. The scaling parameter Uk must be computable on both 
sides of the channel. Since the receiver side knows only jkSk, this parameter must be updated as 

<J k = 5 fc ((7s)o _1 ), <J > Y , 

where (js) k := 7feS fc . 

The transmitted signal Sk is randomly lost due to unreliability in the channel caused by congestion or 
delay. The packet reception/loss state at time k is represented by the random variable jk <S {0, 1}. If jk = 
then the packet is lost at time k; otherwise, it arrives successfully. The process {7fc}^L ^ s independent and 
identically distributed (i.i.d.) with the loss probability specified by p £ [0, 1), i.e., Prob ({7^ = 0}) = p. 

The decoder Dk converts the received signal (js)k to the interval 34 C [— Cfc/2, Ofc/2] byJ4 = Dk ((7s)o,7o)- 
The interval is an estimation set that the plant output y k should be included in. If the packet Sk arrives, 
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i.e., if 7fe = 1, then the output interval 34 corresponds to one of the partitions of the quantizer; otherwise, 
34 is equal to the entire input range of the encoder. Hence, we have 

' [-^ + ft*, + ft (i + 1)) if 7fe = 1 and s fe = i, i £ {0, 1, . . . , TV - 2} , 
^= [t-ff.^] if 7 fe = lands fe = JV-l, (5) 

The controller C fe generates the control input Ufe based on the estimation set 34 as = Cfc(3 7 o )■ 

It is remarked that the scaling parameter should be large enough to cover all possible inputs to the 

encoder. Otherwise, the quantizer may be saturated, in which case we lose track of the plant output yk- 

On the other hand, if we take cr^ large, the quantization error also becomes large. Moreover, to achieve 

stabilization of the system a k should decay gradually. 

We determine the scaling parameter as follows. At time k, the encoder and the decoder predict the next 

plant output Uk+i based on the observed yfi • Let 34+1 C M be the set of all possible outputs yt+i of the 

uncertain system (1). Then the scaling parameter <Jk+i > is chosen such that 

CTfc+i > M34~+i), ( 6 ) 

where fi(-) denotes the Lebesgue measure on M. 

As the prediction set 34~+i, we employ the following one. 

Definition 1 The prediction set of the plant output yt+i constructed at time k is defined as follows: 

34+1 : = WiVk + ■■■ + a'„y'k-n+i ■ a' 1 £A u ...,a' n £ A„, y' h £ 34, ■ ■ ■ , y' k - n +i e 34-n+i}- (7) 
Under this definition, our prediction strategy is to use the iuforma-tiou rega.rd.iug . . . , y&— n+i 

inde- 
pendently such that j/fe-i+i £ 34-i+i for each i = 1,2, ...,n, where 34-i+i is the interval received on the 
decoder side at time k — i + 1. Then, clearly, /ix(347 + i) is large enough to include yt+i, and it is computable 
on both sides of the channel. 

Remark 1 The class of prediction sets in Definition 1 is employed to pursue an analytical approach. There 
may be other prediction methods that generate a less conservative prediction set 34~+i satisfying m(3^jT +1 ) < 
/i(3 7 fc " +1 ). If, for example, we take yk, • • • , yu-n+i not independently but by looking at the correlations among 
them, then the estimation sets 34-i, • • • , 34-n+i from times before k and 347+1 ma y shrink. In the case of 
uncertain plants, however, it is difficult to analytically minimize the Lebesgue measures of the prediction 
sets containing yu+i, though it may be possible numerically [15]. 

The control objective is stabilizing the system depicted in Fig. 1 in a stochastic sense as described below. 

Definition 2 The system depicted in Fig. 1 is mean square stable (MSS) if the plant output yk asymptoti- 
cally goes to zero in the mean square sense for all possible uncertainties within the bounds in (2), i.e., 

lim E[|y fe | 2 ] = 0, Voi, fc £ Ai,Va 2 ,k <^ A 2 , ■ ■ ■ , Va„ ifc e A n , Vfc £ Z+. 

k— >oc 

The problem of the paper is to find limitations on the data rate and the loss probability for the system 
to be MSS. 

3 Necessary Condition 

In this section, we present the first main result of the paper: A necessary condition for the system to be 
MSS. We introduce two notations representing certain lower bounds on the data rate: 

(|A A .| + e„) V^P 



R£) c ■■= log 



v/l-MIAA-l+en) 2 ' 



- iog 2 , ^-^^~p . 

V 1 - P(\^A-\ + e«) 2 - e ny /l-p 
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The following theorem is the necessity result. 



Theorem 1 If the system depicted in Fig. 1 is MSS, then the following inequalities hold: 

R > i?„ cc := max{i?(° c ), (8) 
1-e 2 

P < Pnec := 77T I , ^ J, 0<€„<1. (9) 

(|Aa«| +e«) 2 -el 

One can verify that R ncc is monotonically increasing with respect to |Aa*|, P, and e„ (under the constraints 
of (9)), and similarly p nec is monotonically decreasing. In particular, as expected, more uncertainty in the 
plant will result in higher requirement in communication with a larger data rate and a smaller loss probability. 

A special case of this result is when there is no uncertainty in the plant, in which case we have e„ = 0. 
Then, the bounds in the theorem coincide with the ones given in [8] . 

Proposition 1 ([8]) Suppose that in (3) the exact plant model is known and is time-invariant, represented 
by the triple (A, B, C). Then, the system in Fig. 1 is MSS if and only if the data rate and the loss probability 
satisfy the following inequalities: 



R > R Y (X A ) := log 2 ^f^l , (10) 

V 1 -p\^a\ 2 



1 



pKpyiX^-.^j—. (11) 



Here, denotes the product of the eigenvalues of A, that is, Xa ■— Y\7=i A»(-<4.). 

It is clear that when e„ = 0, then we have that R nec = Ry(Xa*) and p ncc = py(Xa»)- Thus, Theorem 1 
includes Proposition 1 as a special case. For the case e„ > 0, we have the following inequalities 

i? ncc > max R Y {X A ) = R$ c , 

Pncc < min py(Xa)- 

^AeA n 

By these relations, we have that when the plant is uncertain with e n > 0, then even if we assume the most 
conservative plant dynamics, the limitations Ry and py given by Proposition 1 may not satisfy the necessary 
conditions (8) and (9). 

In Theorem 1, the condition e„ < 1 indicates the existence of a threshold on the uncertainty bound, 
over which the system cannot be stabilized with any controller regardless of the size of the data rate in the 
channel. 

Before starting the proof of Theorem 1, we introduce an inequality regarding the set 34+ 1- By the 
definition in (7), the prediction set 34+1 satisfies 

n 

M^k+i) = ^^{Ayk- l+ i) >ii{A n y k - n +\) (12) 

where 

A i y k -i+i ■= {a-t/fe_ i+ i : a- e A, y'k-i+i & 

for i = 1, 2, . . . ,n. Here, the equality holds by the Brunn-Minkowski inequality [16]. 

Proof of Theorem 1: First, we show that the mean square stability of the plant output y k implies that 
the scaling parameter a k is also MSS. For any measurable set 34 C K, it is obvious that 

2 max \y' k \ > (j,(y k ). 
y' k ey k 
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By the definition (5) of the decoder, the right-hand side of the inequality satisfies n{y k ) = o k jN lk . Hence, 
if linife^oo E[|y fe | 2 ] = 0, then lim^oo E[cr^] = holds. 

In the rest of the proof, we study a necessary condition for a k to be MSS. Notice from (6) that a k+ i is 
bounded from below by m(34+i)- Substitution of (12) into (6) yields 

CTfe+l > /J {A n y k -n+l) ■ (13) 

We next evaluate the right-hand side of the above inequality. This term depends on where the intervals 
A n and 34-n+i he with respect to the origin of the real axis. By (2), we have A n — [a* — e„, a* + e n ], but 
(4) implies that it docs not contain the origin. We first consider the case a* > 0. Using basic notions from 
interval arithmetics [17], we obtain 

ii{A n y k -n+i) = a n fi(y k - n +i) + e n /3(y k - n+1 ), (u) 

where /?(•) is given by 

(y k + y k ifo<34, 

y k _- y k iiy^<o<y k , (15) 
{-y k -y k xy k <o, 

y k := sup y' k , y k := inf y' k . 

y' k ey k Vk^k 

Similarly, for the case a* < 0, we have 

M (A n y k -n+i) = -<|u(^- B +i) + e n p(y k - n +i)- (16) 
Clearly, we can write (14) and (16) in one form: 

M {A n y k - n+1 ) = \a* n \n{y k -n+\) + e n /3(y k - n +i)- (17) 
Though the value on the right-hand side of (17) may vary with 34-n+i, the inequality (13) holds for any 

y k -n+i- 

We claim that the maximum of \x (A n y k - n +i) in (17) over all possible intervals 34-n+i at time k — n + 1 
can be written as 

max fi(A n y k - n +i) = VkVk-n+i, (18) 

where r\ k is the random variable given by 

. |<|+max{JV^-"+i-l,l}e n 

% ■= W^Tl ' ( 19 ) 

In (17), the first term is |a* \/-i(y k - n +i) = \a n \ <T k-n+i/N lk - n+1 . For the second term, we must consider the 
following three cases. 

(i) < y k _ n+ \\ In this case, by (5), it is necessary that 7fc-„+i = 1 and N > 2. From (5) and (15), we 
have 

N-l 

max f3(y k - n +i) = — — <T fe _„+i. (20) 

yk-n+l iv 

(ii) y k - n +i < < y k -n+i'- If 7fc-n+i = or N < 2, then this condition is satisfied for any 34-n+i- By 
(15), we obtain (}(y k _ n+1 ) = y k -n+i - 34-n+i = M^k-n+i)- Hence, we have 

max p{y k - n+ i) 
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(iii) 34-n+i < 0: This case can be reduced to (i). We hence obtain (20). 
From the above, the relation in (18) is derived. 
By (13) and (18), it holds that 

o-fc+i > rjkUk-n+1- (21) 

This is a key inequality, which is a consequence of <7k being MSS. Noticing that at > for each time k, take 
the square of both sides of (21) as 

E[al}>n4al_ n+1 }=n4}E[*t_ n+1 }. 

Here, the equality holds due to the independence of <Jk-n+i an d Jk-n+i', note that dk-n+i is an element 
of the cx-field generated by the sequence 7q~". The above inequality shows that if <7fc is MSS then % must 
satisfy E[r^] < 1. 

The next step is to obtain a lower bound on TV from the inequality E[ry^] < 1. For this, we must consider 
two cases. 

(i) 1 < N < 2: By the definition of r]k in (19), in this case we have % = (|o* | + e n )/N lk - n + 1 . Hence, it 
holds that 

„2l /i * I . , \2 



nvi]=P(K\ + e n y + (i- P )(^ ] ^) <i 



(Kl + ^n) 2 " 

Using i? = log 2 iV, we can simplify this as 

i? > 4° e ) c , p < pi°l (22) 
(ii) N > 2: As in (i), E[r^] < 1 is equivalent to 

«>*&■><■&:- K p;4i.. - (23) 

We next simplify these necessary conditions (22) and (23) into one form. First, we consider the constraints 
on the data rate R. According to (22) and (23), a lower bound on R is i? n cc if < R < 1, and R^ec otherwise. 
However, by the definitions of R^cc and i? n cc, we can show the three relations 

R-nec > 1 =^ -Rncc -^noci 

R (o) = R W ^ R (o) = R (i) = l 

-"'nec Jt ncc ^ Jt ncc -"rice i- 

That is, the bound i? ne c is larger than R^Jc when R^cc > 1, is smaller when R^cc < 1, and crosses i? n oc when 
RnJc — I- Thus, the constraints on R can be reduced to the desired bound R > max{i? n el, i? n oc} in (8). 

Next, we turn to the loss probability p. We show that for the case (i), i.e., if < R < 1, then it is clear 
that p < pnoc in (22). Since it holds that R > Rnec, we have i? n ec < 1- Here, let p* denote the p such that 
Rnec = 1 holds. Since R n °ec is monotonically increasing with respect to p, i? n cc < 1 implies that p < p*. 
On the other hand, after some calculation, we can show that p* < pnee- Thus, the condition p < p n °J c is 
automatically satisfied. Moreover, we can also show that p* < p n oc- Hence, p < piec is necessary also in the 
case < R < 1. 

Finally, since pnec must be larger than zero, as a necessary condition, we need < e„ < 1 in (9). ■ 
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4 Sufficient Condition 



In this section, we present a sufficient condition for the existence of a stabilizing feedback control scheme. 
When we consider a practical quantization scheme, we cannot choose the quantization level TV as a noninteger. 
Therefore, we assume TV € Z and TV > 2 in this section. 

Given a certain data rate R, or TV, we employ the control law for the scaling parameter as 



and that for the control input as 



(24) 



(25) 



where ijk-i+x ■= (34-i+i + 34-i+i)/2- Notice that we select the minimum scaling parameter at satisfying 
(6), and take the state feedback control using the nominal values a* and the centers of 34-i+i- 
For i = 1, 2, . . . , n, we introduce the following random variables Oi.k- 



0, 



,fc : 



\a*\+ei 

|<|+ e ,(Af-l) 
N 



if jk-i+i =0, 
if 7fc_ i+ i = l and A^O, 

-{^W^^i} if 7fc-i+i = l and A t 30. 
This can be used to bound the interval A34-i+i as 

Moreover, define the random variable matrix i/r fc containing the random variables . . . , by 

1 ••• 



(26) 





9n,k ^n—l,k 



l 



where Tk := {7fe,7fe-i, • • • ,7fc-n+i}- Here, the process is a Markov chain which has 2" states given by 

r« : ={o,...,o,o}, r< 2 ) :={o,...,o,i}, r< 2 ") : ={i,...,i,i}, 

and the transition probability matrix P e K 2 " x 2 " is given by 



P := 



P 
P 



1-p 
1-p 



P 
P 



1-p 
l-p 



where the element is equal to the transition probability from rW to Furthermore, we define the 
matrix F using Hr k and P: 



where 



F := F1F2 



Fx := P T ® I n 2 , F 2 := diag (H rW H r(1} , . . . , i? r( 2«) ® ff r p«) ) . 



Here, diag(-) denotes a block diagonal matrix and ® is the Kronecker product. The following theorem holds 
by applying results from the Markov jump linear systems theory [18]. 
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Theorem 2 Given the data rate R = log 2 N and the loss probability p£ [0, 1), if 

p(F) < 1, (27) 

then under the control law using (24) and (25), the system depicted in Fig. 1 is MSS, where /?(•) represents 
the spectral radius of a matrix. 

Remark 2 For the special case of scalar plants (n = 1), the condition (27) is equal to the following condi- 
tions: 

R > J?W, p < p noc , < e„ < 1. 

This shows that if the necessary condition in Theorem 1 holds and the data rate satisfies 2 R = N > 2, then 
the sufficient condition is also satisfied. Thus, Theorems 1 and 2 are necessary and sufficient for the case 
n = 1. 



Remark 3 In [12] and [13], sufficient conditions for stabilization of uncertain plants via lossless channels 
(p = 0) are given. For the case n = 1, the condition in [12] is 



R > R P := log 2 



and the one from [13] becomes 



R > R M := log 2 



|AA-|-ei(|A A .|+ei) 
l-e 1 (2|A A .|+2e 1 + l) : 

|Aa-| 



1-ei 



(28) 



(29) 



By Remark 2, our sufficient condition (27) on the data rate equals R > R^Jc for the case n = 1. It is easy 
to verify that J?nec < Rp,Rm- Thus, our sufficient condition is less conservative than (28) and (29). For 
general order plants, it is difficult to compare Theorem 2 with the bounds in [12] and [13] because the types 
of plant uncertainties are different. 

Proof of Theorem 2: We first prove that if E[<r^] — > then E[|j/fc| 2 ] -4- as k — »• oo under the control 
law. By substituting (25) to (1), we obtain 



\yk+i\ 



^ {aiVk-i+i - alVk-i+i) 



- ^2 (l a * _ a *lly*;-i+l| + N \Vk-i+l - Vk-i+l\)- 



(30) 



Regarding the first term in the far right-hand side, we have 

|a» - a* \ < t l (31) 
from the bound on the uncertainty. Moreover, since yt-i+i is the midpoint of 34-i+ii we have 

0~k-i+l 



\Vk-i+i\ < 



(32) 



Similarly, the second term \y/--i+i — yt-i+i\, which corresponds to the quantization error, is bounded by 
(Tk-i+i as 



\yk-i+i - yk-i+i\ < 2fl7fc _, +1+1 < 



(33) 
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Applying these inequalities (31)-(33) to (30), we have 

1 



\vk+i \ < ^2^1 
i=i 



i 

N 



« T \Q>i\ „ 



i=l 



(34) 



where di is defined as follows: 

di := £j 

By squaring both sides of the inequality in (34) and taking expectations, we obtain 



'4 



I + i|ifor» = l,2,...,n. 



E[|y fe+ i| 2 ] <^^^E[a fe 
»=i j=i 



»=i j=i 

Here, the second inequality holds by the fact that er/c > for all k > and the Schwarz inequality. Hence, 
we have that if E[cr|] — > then E[|yfc + i| 2 ] — >• as fc — > oo. 

Next, we prove that the condition (27) implies that a k is MSS under (24) and (25). By (24) and the 
equality in (12), we have 



Cfc+l = ^2 H (A^fc-i+l) • 



(35) 



For the ith term /x (^j^fe-i+i), it holds that 



(|aJ|+ C< )Myfc-i+i) 



if yk-t+i 9 o, 



|aj|/i(y fc _ i+ i) + e\y k - i+ i + y k -i+i\ if ^fc-i+i ? and ^ ^ 0, 



2ei max 



if 34-i+i ^ and Ai 3 0, 



by using basic results in interval arithmetics [17] for i = 1, 2, . . . , n. By taking the maximum of [i (Aiyk-i+\) 
over all possible we have 

M (A^fc-i+l) < ^i.feCfe-i+l, 

where (^fc is given in (26). Thus, from (35), it follows that 

Cfe+l < ^l.feCTfe + #2,fcCTfc_l + • • • + 9 ni kCTk-n+l- 

This can be expressed as 

Cfe+i < -ff/cCfej 

where Cfc : = [cfe-n+i • • • Ofc] T and the inequality is in the element-wise sense. From [18, Theorem 3.9], it 
follows that the inequality (27) implies that £fc and hence a k are MSS. ■ 
We now illustrate the theoretical bounds on the data rate obtained in Theorems 1 and 2 by a numerical 
example. Consider an uncertain plant of second order, where a\ = 1, ei = 0.05, ti = 0.05, and p — 0.05. 
In Fig. 2, we plot the theoretical bounds on the data rate R versus |Aa*| = 1^2 1- The vertical dash-dot line 
represents the supremum of \\a*\ such that p < p ncc holds. Hence, the necessary condition (9) does not hold 
on the right side of this line. The figure shows a certain gap between the necessary condition (solid line) and 
the sufficient condition (dotted line). However, the gap is fairly small when |A,4*| is sufficiently smaller than 
the critical vertical line and is about one bit. 
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Figure 2: Data rate limitations: The necessary condition (solid), the sufficient condition (dotted), and the 
sup \\a»\ s.t. the loss probability p = 0.05 satisfies p < p ncc (dash-dot) 



5 Limitation on the Average Data Rate 

In this section, we introduce a time-sharing protocol, and then study conditions on the average data rate to 
stabilize the system. So far, we have considered that the quantization level N must be taken as an integer. 
Hence, for a data rate R satisfying the sufficient condition in Theorem 2, the minimum quantization level 
is given by N = \2 R ~\, where [•] denotes the ceiling function. Thus, the quantization level may be larger 
than 2 R due to the integer-constraint. In a time-sharing protocol, the observation of the plant output is 
m-periodic, where m G N is the duration of the time cycles. Hence, by employing the protocol we can choose 
the quantization level N m — \2 R ~\ and hence a nonintcger N as the average quantization level. 

It is known that when we know the exact plant model (A, B, C), by taking the duration m large enough, 
we can achieve the necessary bound in Theorem 1. That is, for any data rate greater than the bound, there 
exists a control law with a feasible quantization level [4,5,8]. 

Here, we will see that also when the plant is uncertain, increasing the duration m may help to reduce 
the required data rate. However, for the duration above a certain value, we cannot stabilize the system even 
with infinitely large data rate. This is because a long duration, i.e., low observation frequency, causes the 
accumulation of error in the state estimation due to the uncertainty. 

To simplify the discussion, in this section we consider the scalar plant given by 

Vk+i - akUk = Ufc, (36) 
a k G A := [a* - e, a* + e] , e > 0, \a*\ - e > 1. 

Here, we introduce the encoding and decoding rules based on a time-sharing protocol. The basic idea is as 
follows: Divide the time into cycles of a constant duration m G N as {mj, mj + 1, . . . , m(j + 1) — 1}, where 
j G Z+. The encoder observes the output y m j at time mj, and sends it to the decoder during the cycle from 
time mj to m(j + 1) — 1. Over the lossy channel, the decoder receives RYmLo lmj+i bits in total in this 
cycle, and compute the estimation set y m j of y m j at time m(j + 1) — 1. 

More specifically, the encoder and the decoder function as follows. At time mj, the encoder quantizes 
y m j with the quantization level N m . We use the uniform quantizer as in Section 2. The quantized signal 
is divided into packets of R bits, and is then transmitted to the decoder. If a packet drops out then the 
encoder retries to send the same packet. Notice that the encoder has access to the loss states because of 
acknowledgement signals from the decoder. 

The scaling parameter <ik is taken to be large enough to cover all possible outputs at the next step. In 
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this section, since the observation is m-periodic, we slightly change (6) to 

v m (j+i) > V(y~(j +1) ), (37) 
where is the one cycle prediction set based on the estimation set y m j, and is defined as follows: 

y~ [j+1) ■= {{aTy'mj + (aT^Umj + (a') m - 2 u mj+1 + ■ ■ ■ + a'u m{j+1) _ 2 :a'€A, y' mj e y mj } • (38) 

Note that we have to take account of the control inputs Uk from k = mj to m(j + 1) — 2 since they affect 
the length of the prediction set. 

We now study a condition for the stabilization of the system under the protocol described above. 

As we have shown in the proof of Theorem 1, for lim^oo E[|j/fc| 2 ] = 0, it is necessary that lim^oo E[ct 2 ] = 
0. To derive a condition for stabilizing o^, we follow an approach similar to that in Section 3. In view of 
(37), we evaluate the length of the prediction set D^CH-i)' By (38), we have 

My-y+ij) > M ({(aTy'mj ■■ a' € A, y' mj e y mj }) (39) 

= Ua*\ m ii(y mi ) + 5 + \y mj \ + 5-\y mj \ ity mj Jo, 

\(\a*\ m + 6+)fi(y mj ) if 3^ 9 0, 

where 5+ and <5_ are nonnegative values defined as <5 + : = ( | a* | + e) m — \a*\ m and 6- := — ( | a* | — e) m + | a* | m , 
respectively. The maximum of the right-hand side of the above equality with respect to y m j can be expressed 
as 

f ]a *l7 +6+ °-mi if 1 < M mi < 2, 

max/z {{{aTv'mi ■ a' ^ A, y' mj € y mj }) = , ^_ s _ s++s _ s (40) 
y m] y y-—— + ___ j Umi if M m] > 2, 

where M m j is given by 

M mj := N^^^+K 

Notice that the right-hand side of (40) can be written as K m j<j m j, where the n m j is the random variable 
defined as 

^" := ^( |a * r + max {^' 1 } ,5+ + max {^T i " 1 ' }' 5 - 

Thus, by (37), (39), and (40) it follows that 

(41) 

Since <j m j > for all j, and n m j is independent of a m j, a k is MSS only if R := E[/c^,-] < 1. Thus, we 
have arrived at a necessary condition for the overall system to be MSS. Note that the expectation R is 
time-invariant because it depends only on the process {7fc}^ , which is i.i.d. 
As a stabilizing control law, we employ the following: 

o-mtf+i) =n(y~ u+1) ), (42) 

Jo if*e{0,l,...,ro-2}, 
u m j+i = < . Mm „ . (43) 

[-{a*) m y mj if z = m — 1, 

where y mj := (y mj +y mj )/2. 

The proposition below presents a necessary and sufficient condition for the stabilization of the system 
employing the time-sharing protocol. 



12 




J \ 2 3 4 

Duration of cycles 

Figure 3: Limitations on the average quantization levels: The necessary bound (solid), supm s.t. (6+ + 
S-)/2 < 1 holds (dash-dot), and the sufficient and feasible bound (cross) 



Proposition 2 Consider the system depicted in Fig. 1 with the uncertain scalar plant (36), and a control 
law under the time-sharing protocol. If the system is MSS, then the following inequality must hold: 

k < 1. (44) 

Conversely, if (44) holds then the system is MSS under the control law (42) and (43). 

Proof: (Necessity) Refer to the discussion before the proposition. 

(Sufficiency) As in the proof of Theorem 2, we have that if lim^oo E[af\ = then limfc_ ) . 00 E[|yfe| 2 ] = 
under the control law (42) and (43). Moreover, equality holds in (39) and hence in (41). Thus, for any N 
and m that satisfy (44), N m <G Z, and N m > 2, we have that a k is MSS. ■ 

From Proposition 2, we can obtain the explicit limitation on the data rate when the channel is lossless, 
i.e., for the case p = 0. To describe the limitation, let 

— 1 la* l m — S 

fl - : ^m l0g2 l- (( 5 + + M/2 - 

Then, the condition (44) is equivalent to 

R > R DCC := max {rW, } , < < 1. (45) 

Here, i?noc is the necessary bound defined in Section 3. Clearly, the bound on the average data rate depends 
on the duration m, and is equal to i? ncc if m = 1. In particular, we cannot stabilize the system even with 
large R when m is greater than a certain value; if m is large enough to satisfy (S + +S-)/2 > 1 then R^ec = oo- 
We illustrate the limitation on the average quantization level given in Proposition 2. Consider an un- 
certain scalar plant, where a* — 3.3 and e = 0.025 with a lossless channel p = 0. In Fig. 3, we plot the 
limitations on the average quantization level N versus the duration m of the cycles. The vertical dash-dot 
line represents the supremum of m such that the condition (S + + <5-)/2 < 1 in (45) holds. The necessary 
bound (solid line) increases with respect to the duration and becomes infinitely large as it reaches the dash- 
dot line. The cross marks represent the sufficient and feasible average quantization levels, i.e., the levels 
satisfying (45), N m £ Z, and N m > 2. We can find that the necessary bound on iV takes its minimum at 
m = 1. However, as the cross marks show, when we take account of the feasibility the duration m = 2 is the 
best. 
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6 Conclusion 



In this paper, we have studied the stabilization problem of uncertain networked control systems under the 
presence of data rate constraints and packet losses. In particular, we have derived a necessary condition and 
a sufficient condition for the stability of the closed-loop system. These conditions highlight limitations on 
the data rate, the loss probability, and the uncertainty bounds, and generalize existing results on nominal 
plants to the case with parametric uncertainties. 
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